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Abstract 

We investigate the effects of homogeneous general dark energy on the non-linear mat¬ 
ter perturbation in fully general relativistic context. The equation for the density contrast 
contains even at linear order new contributions which are non-zero for general dark en¬ 
ergy. Taking into account the next-leading-order corrections, we derive the total power 
spectrum in real and redshift spaces. We find that the observable galaxy power spectrum 
deviates from the ACDM spectrum, which is nearly identical to that in the Einstein-de 
Sitter universe, and the relative difference is about 10% on a scale of the baryon acoustic 
oscillations. 



1 Introduction 


After the discovery of the accelerated expansion of the universe from the observations on distant 
type la supernovae P, among the biggest puzzles in cosmology is the identity of dark energy, the 
“driving force” of the acceleration that occupies about 68% of the energy budget of the present 
universe [2J. The simplest possibility is a small but non-vanishing cosmological constant [3], 
whose right amount would be then explained possibly by some unknown symmetry of quantum 
gravity or anthropic principle based on the landscape of vacua. Dynamical alternatives are 
non-trivial scalar fields such as Galileon [4], or modification of gravity like massive gravity [5] 
or f(R) theories [6j. 

While the identity of dark energy has a variety of theoretical origins, its properties are 
constrained by observations [2], from which we understand how dark energy has influenced the 
history of the universe: at very early times, for example when the cosmic microwave back¬ 
ground (CMB) was generated, dark energy was completely negligible. At later stage when 
non-linearities in cosmic structure are developed, dark energy becomes significant and affects 
the evolution of gravitational instability. Therefore, the effects of dark energy should emerge 
more prominently at non-linear level. Possible signatures of dark energy in non-linear cosmic 
structure can be measured with more accurate data from large scale galaxy surveys into deeper 
redshift in near future such as HETDEX [7], DESI [B], LSST |Tj and Euclid [1U], which are 
to probe very large volume comparable to the horizon scale. This also demands fully general 
relativistic approaches beyond the Newtonian theory. 

In this article, we study how homogeneous dark energy with a time-varying equation of 
state modifies non-linear matter power spectrum in the comoving gauge. The same topic was 
studied in the Newtonian perturbation theory but the effects of general dark energy are totally 
insignificant HU; see also |l2j . We find that in the presence of dynamical dark energy the 
next-to-leading corrections to the power spectrum can give rise to a deviation as large as a 
few percent from both the non-linear power spectrum in ACDM and the counterpart in the 
Einstein-de Sitter (EdS) universe at scales close to the baryon acoustic oscillations (BAO); 
meanwhile the power spectrum in ACDM is almost identical to that in the EdS universe with 
less than 0.1% difference. This notable difference between the power spectrum in fully general 
relativistic context and the corresponding counterpart in the EdS universe, which is usually 
taken in the Newtonian theory, is thus a sharp feature of general dark energy models. For the 
power spectrum in the redshift space which is directly observable, the deviation becomes more 
prominent. At BAO scales the deviation from ACDM reaches about 10% or so. 

We also investigate how the homogenous dark energy affects the pure relativistic contri¬ 
butions coming from the curvature perturbation, which is heavily suppressed in the EdS uni¬ 
verse mm- The curvature effects decrease more than 100% from the EdS case in the presence 
of dark energy, including a cosmological constant, while their contribution is not significant at 
all. 


2 Modified equations with general dark energy 

We consider a flat Friedmann-Robertson-Walker universe as our background. The energy den¬ 
sity of the universe consists of pressureless matter and non-interacting dark energy, which is 
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a good approximation for the epoch of the evolution of cosmic structure after the generation 
of CMB until now. Furthermore, to simplify the analysis, we assume that dark energy is so 
homogeneously distributed that its perturbation is negligible. That is, dark energy contributes 
only to background; see however [H>] for possible significance of dark energy perturbation. To 
fix the coordinate system, we choose the comoving gauge where T°j = 0 with the spatial metric 
gij = a 2 (r])(l + 2 ip)8ij. Here, dr / = dt/a is the conformal time. Note that the governing equa¬ 
tions in the comoving gauge is known to coincide exactly with the Newtonian hydrodynamic 
equations up to second order in a zero pressure fluid [[[6]. This is, however, not the case when 
we consider dark energy, because of the pressure associated with dark energy which plays a 
genuine relativistic role even at background level D21- For example, the energy conservation 
equation has an additional contribution 

- Kp m ( 1 - A) with A = (1 + w) ^1 <- . (1) 

Here, k is the perturbation in the extrinsic curvature, an overbar means background, and w is 
the equation of state of dark energy. A is thus an indicator of the effects of dark energy and is 
a key parameter in this work. Note that it vanishes in both EdS universe and ACDM. 

The equation for the matter density contrast 5 = 8p m Jp m is then obtained by combining 
the energy conservation and Raychauclhuri equations as 

+ (n + s' -1( 1 , - a )u 2 n m 5 = j\f N +j\f tp +My + a r x , ( 2 ) 


where a prime denotes a derivative with respect to the conformal time and H = a'/a. Note that 
on the left hand side, where we have placed linear terms, we find new contributions with A that 
are absent in the conventional Newtonian perturbation theory. The terms on the right hand 
side are non-linear sources: A/jv denotes the non-linear Newtonian source Em, and the pure 
relativistic effect M v originates from the curvature perturbation ip [EUE]. The remaining two 
terms, A fy and A/a, represent respectively the evolution of p and the departure from ACDM 
and EdS, and are given by 
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where (3 r = A -2 (7 Af3Ap/2 + 2 P ,l Ap^ — /3’ y V,p/2 + 2<pA 2 /3 + 3<p ,l AP ti ). Here, P is the scalar 
perturbation in the 0 i component of the metric given by g 0 j = a 2 (p)/3j, A = 8 l Widj is the 
spatial Laplacian and A -1 is the inverse Laplacian operator. It is important to note that while 
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the effects of the curvature perturbation appear only from third order US). general dark energy 
modifies the equation from linear order as can be seen from (J2]) and Q. This is because dark 
energy changes the background expansion rate and thus its effects permeate throughout all 
order in perturbation. 

At linear order, the solution of the curvature perturbation ip can be found as 


ip = —A 


-i 


+ ^n 2 n m 5 


n 2 f 

1-A 


1 + |(1 - A )^ 


A <5, 


(5) 


where / = dlog-Di/dloga and Di is the linear growth function. In the EdS universe, the 
expansion rate of universe is same as the inhomogeneity growth rate at linear order: Di = a 
so that T-[ 2 Di is constant. Therefore, ip is time-independent in the EdS universe. Similar 
conclusion of a constant ip can be drawn in ACDM. On the other hands, ip decreases when 
we consider general dark energy because of the faster expansion rate and slower growth of 
inhomogeneity. Thus, ip' also indicates the effects of general dark energy. 


3 Non-linear power spectrum in real space 

To compute the power spectrum perturbatively, we expand the density and velocity divergence 

perturbations as 5 = <5i + 8 2 + £3 H-and k = + n 2 + k 3 H— •. The second and third order 

solutions are found in terms of the n-th order symmetric kernels F n along with time-dependent 
coefficients as 

<5 2 (fc, a) = D\(a ) Y c 2i {a) J ^^h (3) ( fc ~ 912 ) F 2i (q ± , q 2 )Si(qi)S 1 (q 2 ) , 

5 3 (fc, a) = D\{a) Y c 3 i(a) J ~ ^ ^ ( fc “ F 3*(9 i, <?2, q3)Si(qi)Si(q 2 )S 1 (q 3 ) 

+ Dl(a)n 2 (a) <%i(a) J d ^ _ ^ 123 ) F^( qi ,q 2 , q 3 )5 1 (q 1 )5 1 (q 2 )5 l (q 3 ), 

(6) 

where c ni = D ni /D ”, c^ = D^J (Dfl-L 2 ) and q\ 2 ... n = Qi- The superscript ip means orig¬ 
inated from ip. The solutions for K n /l-l have exactly the same form as S n , while the time- 
dependent coefficients are replaced by k ni = K ni / (HD™) with K ni being Ath growth function of 
K n . In the EdS universe, the coefficients are hxed as certain numbers: c 2a = 3/7, k 2a = -1/7, 
c 2b = 2/7, k 2b = 4/7, c 3a = k 3a = c 3b = k 3b = 1/2, c 3c = k 3c = 2/7, c 3d = -2/63, k 3d = -1/21, 
c 3e = -1/9, k 3e = 8/21, c 3/ = 8/63, k 3f = -2/21, cf a = kf ia = -5/4, cf 6 = -5/14, and 
k^ b = —20/7. Since the universe is close to EdS when 6 begins to grow, these numbers are set 
to be the initial conditions for these coefficients and we evolve them. Note that in S 3 , unlike 
ip, we cannot separate the contributions originated from A because as stated in the previous 
section dark energy changes the background expansion rate. Also note that the contribution 
from ip carries H = aH, which can be identified as the comoving horizon scale kn- 

The non-linear power spectrum is defined by ( 8(k)8(k')) = (2 tt) 3 5^ 3 \k + k')P(k). Under 
the perturbative expansion of 8, we can write P = Pu + P 22 + P i3 + • • •, where the subscript 
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denotes the order of perturbations being correlated, so that Pn is the linear power spectrum 
between two hi’s and P 22 + -P 13 is the next-to-leading order, viz. one-loop power spectrum. The 
power spectrum with one-loop corrections is 


P(k,a) = D 2 1 (a)P 11 (k) 
k 3 
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(7) 


Here, r and £ represent respectively the magnitude of the dummy integration momentum by 
q = rk with 0 < r < 00 and the cosine between k and q by k- q = kqx = k 2 rx with —1 < x < 1. 
Note that in the EdS universe and ACDM, P 33 becomes purely Newtonian and Pf ?j alone carries 
the relativistic effects ra 

To evaluate the power spectrum ([7]) for general dark energy, we parametrize the equation 
of state as [19] 

w(a) = w 0 + (1 — a)w a . (8) 

Then we calculate the linear power spectrum P n from CAMB [2D] with the best-fit parameters 
from the Planck 2015 data [2], which is used to evaluate the next-to-leading contributions 
P 22 + P 13 + P \3 in (J7|). In Figure 0 we present in the top panels the total power spectra for 
different dark energy models at z = 0 along with the one-loop components P 22 + Pi 3 and Pf s - 
In the left (right) column, with a fixed value of w a = 0 (w 0 = —1) we vary w 0 (w a ): w 0 = —0.8, 
—1.0 (ACDM) and —1.2 with w a = 0 in the left panel, and w a = 0.5, 0 (ACDM), —0.5 and 
—1.0 with wq = —1.0 in the right panel respectively. In the bottom panels of Figure [l] the 
relative differences of the dark energy models from ACDM are presented. The deviation from 
ACDM is not significant and nearly constant on large scales (k < O.lh/Mpc), but it increases 
from k ~ 0.1/i/Mpc and the relative difference becomes more or less 10% at k ~ 0.2/i/Mpc and 
even larger on smaller scales. We can note that the power spectrum for Wq > — 1 (wo < — 1) in 
the left panel and w a > 0 ( w a < 0) in the right panel is smaller (larger) than that in ACDM. 
Since we normalize D m o as the Planck 2015 best-fit value, the energy density of dark energy 
increases (decreases) when w < — 1 (w > — 1). Thus the matter-dark energy equality is later 
for w < — 1 than for w > — 1, although the expansion rate for w < — 1 is faster than that for 
w > — 1 at present. As a result, the inhomogeneity can grow more for w < — 1 than w > — 1. 
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Usually, one uses the non-linear power spectrum in the EdS universe to an arbitrary cosmo¬ 
logical model by replacing the scale factor with the linear growth factor of that model, while 
leaving the time-independent part of the spectrum identical to the EdS one: 

P(k,a) = Dl(a)P u (k) + D\(a) \P 22 (k) + P 13 (*;)] EdS . (9) 

Compared with the true power spectrum (J7]), this is obviously inconsistent because in general 
dark energy models the non-linear kernels Fi s as well as the time-dependent coefficients q’s of 
the non-linear solutions (Jb]) are different from the EdS counterparts. In Figure [2j we compare 
the power spectra for different dark energy models with the corresponding counter EdS power 
spectra (§. As in Figure[lJ in the left (right) panel we fix w a = 0 (wq = —1) and vary Wq (w a ). 
The top panels show the total power spectra, while in the bottom panel the relative differences 
of the dark energy models from the corresponding EdS counterparts are presented. Since both 
J\f^' and A/a vanish for ACDM, the relative difference of the total power spectrum in ACDM 
from the corresponding EdS one is only from Af v and is very small, O (10~ 2 — 10 _1 )% over all 
scales [13]. On the other hand, the power spectrum in general dark energy differs substantially 
from the EdS counterpart because of the large Af\ contribution: the relative difference notably 
increases from k « O.lh/Mpc and becomes larger than 10% for k > 0.2/i/Mpc. While even 
larger deviations on smaller scales may not be reliable because of strong non-linear effects, the 
analytic non-linear power spectrum (J7|) gives very good agreement with numerical simulations 
in weakly non-linear regime [2T]. 

It is also interesting to take a closer look at P'y A . In the top panels of Figure [lj we separately 
show Pf A for different dark energy models. While P 22 + P 13 contributes more significantly on 
smaller scales k > O.lh/Mpc, Pf 3 extends to larger scales. Furthermore, even for ACDM the 
relative difference from the counter EdS power spectrum is larger than 100%, i.e. always smaller 
in the presence of any form of dark energy on all scales. Thus, Pf A is strongly influenced by 
dark energy and is in principle interesting, but its magnitude is heavily suppressed. 

4 Non-linear power spectrum in redshift space 

The power spectrum presented in the previous section is computed in the real space. In reality, 
however, we identify the position in terms of redshift. Thus to extract observational prospect, 
we need to consider the power spectrum in the redshift space which is directly related to 
observations. 

The proper transformation from the real space to the redshift space includes all relativistic 
effects such as gravitational lensing and the Sachs-Wolfe effects. Note that although the real 
space solutions in the previous section are not gauge-invariant, those in the redshift space 
obtained from this transformation are gauge-invariant 12 a. In this article, we take into account 
only the Doppler effect, because the other relativistic effects are suppressed at the BAO scales 
in which we are interested [221123]. Then the density contrast in the redshift space is written 
as, up to third order [24]. 

6 S = 5 r - <9 ||U + (<9 ||Uf - 5 r d\\U - A Ud\\5 r + A Ud\U + 5 r {d\\U) 2 

- (d\\Uf - 3AD<9||/7<9jU + 2/\Ud\\Ud\\8 r + A U5 r d^U + ^(AU) 2 d^S r - ^(A[/) 2 <9fZ7, 

(10) 
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where S s (S r ) means the galaxy density contrast in the redshift (real) space. With a constant 
bias factor b, we can write 5 r = b5 m with S given by ([6]) . Since in the comoving gauge 
the coordinate time is equivalent to proper time, the bias factor may show up as it is without 
further manipulations [26] • The potential U (r) is defined by U = h ■ v/Ji, where v and h are 
respectively the peculiar velocity and the unit vector along the line-of-sight direction (LoS). 
Since k = —V v, U can be written as U = — d\\ where d\\ = n-V is a partial differential 

operator along LoS. 

Then, we can write the observable galaxy power spectrum in the redshift space as 

P»(k, Hi a) = P sll (k, Hi a) + P S 22 (k, h, a) + P s 13 (k, /j, a ), (11) 

where /i = h ■ k/k is the cosine between h and k. Note that GD is no longer isotropic even 
though the power spectrum in the real space is so. The power spectrum along LoS (p = 1) is 
dominant, while the one with /i = 0, i.e. perpendicular to LoS is the same as ([7]) . What is 
important is that hence the deviation from ACDM becomes more significant for the LoS power 
spectrum. In Table [lj we give the relative differences of various dynamical dark energy model 
from ACDM. As can be seen, the deviation is even more enhanced compared to the real space 
power spectrum discussed in the previous section, as large as 10% at around BAO scales, well 
within the expected observational accuracy for future surveys [27]. 


w a = 0 and varying wq 

k [h/ Mpc] 

wo = —1.2 

wo = —0.8 

0.1 

6.8% 

-10.2% 

0.2 

11.6% 

-15.0% 

0.3 

16.0% 

-19.4% 


Wq = — 1 and varying w a 

k [h/ Mpc] 

W a = -1.0 

Wa = -0.5 

W a = 0.5 

0.1 

9.5% 

5.8% 

-11.5% 

0.2 

14.9% 

8.8% 

-15.3% 

0.3 

20.1% 

11.6% 

-19.0% 


Table 1: Relative differences between various dark energy models and ACDM P s /Pacdm — 1[%] 
with /i = 1, z = 0 and b = 2. 


5 Conclusions 

In this article, we have studied the non-linear power spectrum in the presence of general, 
homogeneous dark energy in the context of fully general relativistic approach. We have derived 
the power spectra in the real and in the redshift space. We have found that in general dark 
energy different from a cosmological constant results in substantially large deviations from the 
conventional expectation in the observable power spectrum, as large as 10%, via the next-to- 
leading corrections on scales k > 0.1/r/Mpc.This is quite close to the BAO scales, thus our 
finding should be useful to constrain the properties of dark energy with the large scale BAO 
surveys. 
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